Radio- frequency discharges in Oxygen. Part 1: Modeling 
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In this series of three papers we present results from a combined experimental and theoretical 
effort to quantitatively describe capacitively coupled radio-frequency discharges in oxygen. The 
particle-in-cell Monte-Carlo model on which the theoretical description is based will be described 
in the present paper. It treats space charge fields and transport processes on an equal footing 
with the most important plasma-chemical reactions. For given external voltage and pressure, the 
model determines the electric potential within the discharge and the distribution functions for 
electrons, negatively charged atomic oxygen, and positively charged molecular oxygen. Previously 
used scattering and reaction cross section data are critically assessed and in some cases modified. 
To validate our model, we compare the densities in the bulk of the discharge with experimental data 
and find good agreement, indicating that essential aspects of an oxygen discharge are captured. 

PACS numbers: 52., 52.35.Tc, 82.33.Xj, 47.ll.-j 



I. INTRODUCTION 



Electro-negative gases, such as SFg, CI2, CF4, and O2, 
play an important role in plasma-assisted materials pro- 
cessing. Most notably O2, either as the main feedstock 
gas or as an admixture to halogen- or silicon-based gases, 
is of vital importance for a large variety of etching and 
thin-film deposition techniques [H, [3, % i, i, 0, Q. The 
requirements on the controllability and predictability of 
these processes are so high that further advancement 
of this technology will depend on modeling tools which 
go beyond the macroscopic, fluid-type approximations, 
which, for instance, cannot reliably predict the velocity 
distributions of the species. 

Although electro-negative gas discharges have been 
studied for a long time [S, Ailfllj with significant progress 



made during the last two decades 
plete quantitative description is still lacking, in partic- 
ular, with respect to discharge profiles, structuring, and 
operation regimes l^, lB| . This is not surprising because 
the occurrence of negative ions leads to abrupt changes in 
the ion density (density fronts [17]) which in most cases 
force the discharge to stratify into a central quasi-neutral 
ion-ion plasma and a peripheral electro-positive electron- 
ion plasma [13, [HI, [3, [ll Q El, [H . The transition 
between the two is rather subtle. It can be, for instance, 
accompanied by a double layer (internal sheath [29|). 
Electro-negative gas discharges are thus rather complex 
and the investigation of the spatio-temporal structure of 
the discharge as a function of external control param- 
eters (current, voltage, frequency, pressure, and geom- 
etry) is a g reat theoretical 12 iL \2% l23l . |24 1 and experi- 
mental [11 m, [13, m, [H, [3^ challenge. In addition, 
electro-negative processing gases are reactive molecular 
gases, with internal degrees of freedom, which lead to a 
quite involved plasma chemistry. 



Similar to the kinetic description of a mixture of re- 
acting gases [31|, the modeling of a gas discharge could 



be based on a coupled set of Boltzmann equations for 
the distribution functions of the species which need to 
be treated kinetically. This set has to be augmented by 
Maxwell's equations, or parts of it, depending on how 
the discharge is electrically driven. Even for simple re- 
active gas discharges, with only two negatively and one 
positively charged species, this approach is not practi- 
cal. More promising are methods which track repre- 
sentative samples of simulated particles subject to (non- 
reactive and reactive) collisions as well as electromagnetic 
fields, which are again determined from the relevant parts 
of Maxwell's equations. Although these approaches are 
closely related to the direct simulation method success- 
fully employed for the modeling of the flow of rarefied 



imploy 

gases 3j, |33|, , in the context of gas discharges [35| , it 



is more common to refer to them as particle-in-cell Monte 
Carlo collision (PIC-MCC) methods [H, [13, H HI , be- 
cause without collisions the methods collapse to the PIC 
approach 

, S3] for the solution of the Vlasov problem. 

This is the first paper in a series of three where we re- 
port the results of a combined experimental and theoret- 
ical study of the interplay between plasma-chemistry and 
electrodynamics in capacitively coupled radio-frequency 
(rf) discharges in oxygen. Oxygen is a weakly electro- 
negative gas whose plasma chemistry is strongly affected 
by vibrational, rotational and meta-stable states. Up to 
75 reaction and scattering processes have been listed to 
potentially affect the properties of the discharge ^3] ■ Not 
all of them can be equally important. The challenge is 
therefore to identify (via comparison between experiment 
and modeling) the subset of collisions responsible for the 
experimental findings under consideration. 

In the present paper, we describe such a reduced 
model. It is tailor-made for the investigation of the 
charged species of the discharge, including the formation 
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FIG. 1: Schematic geometry of the rf discharges used in 
Refs. [27l . [30| . We use a planar, one-dimensional model to 
simulate the central axial part of the discharge (thick solid 
line) . 



of ion density fronts, internal sheaths, and the resulting 
stratification of the discharge. We critically assess the 
cross sections used to characterize the selected collision 
channels and point out some inconsistencies in previously 
used cross section data. The validity of our model is ver- 
ified by a comparison of calculated bulk densities with 
experimentally measured ones j27| . The following two 
papers discuss, from the experimental [s^l and theoreti- 
cal point of view, respectively, the sheath region of 
the discharge focusing on the 844 nm double emission 
layer in front of the powered electrode and a comparison 
of experimental and simulated distribution functions for 
positive ions. 

The next section describes the main features of our 
PIC-MCC implementation 



44, m 



It is an exten- 
sion of a code, which was originally designed for the 
particle-based modeling of sheaths in magnetized plas- 
mas 46|, |47| , to the case of electrically driven discharges 
in reactive molecular gases. We focus, in particular, 
on the handling of collision processes, which deviates 
from the treatments usually employed in the plasma 
context [H, 38, [s^ and is closer to the simulation ap- 
proaches used in rarefied gas dynamics 32, Isj]. The 
collisions included in our model are discussed in Section 
III. Where necessary, we combine measured cross sec- 
tions [islligl, [sol, [sli [saJHiJsiJHsl, B [HsUHqI, [eoj with 

model cross sections (6ll |62|. |64 1651 1 to characterize 
collisions and fix the free parameters through a compari- 
son with experimental data. In Section IV we show that 
once the parameters are fixed, good agreement between 
simulation and experiment can be achieved in the bulk 
of the discharge. Section V summarizes the essentials of 
our model and concludes with a short outlook. 



II. METHOD OF SIMULATION 

We are only interested in the central axial part of the 
rf discharges described in Refs. 27, [sO]. Ignoring the 



electric asymmetry between the powered and grounded 



electrode, this part of the discharge can be simulated by 
the planar, one-dimensional model shown in figure[TJ The 
planar model retains only one spatial coordinate x, with 
< a; < L, where L is the distance between the elec- 
trodes, but it keeps the full three-dimensional velocity 
space for the particles (IdSv model). One of the elec- 
trodes is driven by the rf voltage U{t) = Urf sin(27r/r/t), 
with frf — 13.6 MHz and Urf ranging from 75 V to 
800 V, while the other is grounded. Both electrodes 
are assumed to be totally absorbing; secondary electron 
emission is neglected. To mimic the constant oxygen flow 
through the discharge chamber, we enforce in the simu- 
lation volume a constant oxygen pressure ranging from 
10 Pa to 100 Pa. 

The simulation volume of the planar model is of course 
not well defined because the lateral "cross-section" A is 
a free parameter. This parameter actually controls the 
weight of the simulated particles, that is, the number of 
real particles represented by one simulated particle. We 
arbitrarily choose A = A|)g with Xjje the electron De- 
bye length of a reference electron system (RES) which 
we use to initially set-up the length and time scales of 
the simulation (see below). The RES is specified by a 
density n^es and a temperature Tres- Both have to be 
adjusted to the particular experimental conditions. For 
the experiments 27, Is^l, rires ~ 10* - 10^ cm"'^ and 
ksTj-es ~ 10 eV, resulting in approximately 10^ simu- 
lated particles. 

The complex plasma-chemistry of oxygen gives rise to a 
large number of non-reactive and reactive collisions [i^ . 
In tableUwe show the ones with the largest cross sections. 
They are included in our model and will be discussed in 
detail in Section III. We treat only three species kinet- 
ically: electrons (e), negatively charged oxygen atoms 
(0~), and positively charged oxygen molecules (Oj). 
Neutral particles appearing either as educts or products 
in tableware not explicitly simulated. Molecular oxygen 
in its ground state, O2 (feedstock gas), is modelled as 
a reservoir characterized by a pressure p and a tempera- 
ture T. Whereas atomic oxygen, O, vibrationally excited 
oxygen molecules, 02(1'), the Rydberg state of the molec- 
ular oxygen, 02(Ryd), and the molecular meta-stables, 
02{a^Ag) and 02(a^Sg), are only indirectly accounted 
for in as far as their production results in an energy loss 
for electrons. 

The meta-stable 02(0^ Ag) requires special attention 
because it also appears as an educt in the entrance chan- 
nel for associative detachment (17). Its concentration is 
therefore important, and we should actually build-up the 
02(0^ Ag) distribution function, that is, we should also 
simulate 02(0^ Ag) particles. In that case, however, not 
only their volume production process (11) but also their 
volume and surface loss processes should be included. 
This is beyond the one-dimensional model. To take asso- 
ciative detachment, which is known to be an important 
process [23,|6l,[67, 68 1, nevertheless into account, we use 
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FIG. 2: Graphical representation of our implementation of 
the PIC-MCC scheme. 



instead a simple model with one free parameter, which 
can be interpreted as the 02(a^Ag) to O2 density ratio 
(see below). 

The three-species PIC-MCC model describes the 
physics of an O2 discharge by tracing the spatio-temporal 
evolution of an ensemble of charged particles (e, 0~, and 
O^). As usual, to track the velocities and positions of 
the simulated particles, the simulation domain < x < L 
and the time are discretized and free flights are decou- 
pled from collisions (see figure [2]). In order to have stable 
free flights of charged particles (in the mean electric field 
given by the Poisson equation), the spatial resolution Ax 
has to be less then the electronic Debye length, A oe , and 
the time step At should resolve the electron plasma os- 



cillation whose frequency is ujpe [40|, |41| . Both Xjje and 
LUpe are set by the RES. 

The time interval over which collisions and free flights 
are decoupled is of the order of the mean free time, which, 
in general, is much larger then the fundamental time step 
At enforced by the electric field. We distinguish two 
groups of collisions, each of which characterized by an 
average mean free time. The first group of collisions, oc- 
curring on a time scale Atc,i comprises all collisions of 
table [J except Coulomb collisions between charged par- 
ticles, which take place on a different time scale Aic,2- 

To build up a representative sample for the charged 
particles, we first sample the initial distribution func- 
tions for e , ~, and and then repeat the following 
procedure [45| (see figure [2]): 

(i) For m — Atc^i/At time steps, the coordinates and 
velocities of all simulated particles are moved ac- 
cording to Newton's equations of motion. The 
force acting on particles with charge qi is given 
by Fi = —QiV^ where $ is the electric potential 
satisfying the Poisson equation with the boundary 
conditions specified above. Force interpolation and 



charge assignment are done within a linear interpo- 
lation scheme. 

(ii) After ni time steps, a representative sample of 
the first class of collisions is executed, with post- 
collision velocities stochastically determined in ac- 
cordance with energy and momentum conservation, 
followed by free flights with ni — 71,2 time steps, 
where 712 = Aic,2/At, after which a representative 
random sample of Coulomb collisions is performed, 
again consistent with energy and momentum con- 
servation. 

After typically 10^ rf cycles quasi-stationarity with re- 
spect to the rf cycle is achieved, that is, within the sta- 
tistical noise, macroscopic quantities, for instance, den- 
sity and potential profiles, do not change anymore when 
averaged over a rf cycle. 

Our treatment of collisions utilizes a number of con- 
cepts originally developed for the direct simulation Monte 
Carlo of rarefied gas flows [sj 



. Since it differs 



from the handling of collisions in the standard PIC-MCC 
schemes, which is based on the null-collision method [s^ . 
we give some more details. After a free flight is com- 
pleted, we sort the simulated particles into cells (see fig- 
ure [2]). Simulated and feedstock gas particles located in 
one and the same cell (collisions are local) have then the 
opportunity to sequentially perform all possible types of 
collisions with each other. Simulated particles produced 
in the exit channel of a collision are readmitted for the 
next collision in the list. The ordering of particles into 
cells has to be repeated after each collision producing 
or annihilating simulated particles. Thus, as in the real 
system, simulated particles may perform more than one 
type of collision within the time interval over which the 
simulation decouples free flights and collisions. 

The probability of a simulated particle i to make a 
collision of type j is given by 



(1) 



It depends on the density nj of the collision partners 
which act as the targets, the total cross section aj for 
the collision, the relative velocity Uj between particle 
i and a particular collision partner, and the mean free 
time Tj . As can be seen in table HI except for associa- 
tive detachment (17), which we discuss in more detail in 
the next section, the three species model recruits colli- 
sion partners either from the simulated particles or from 
the feedstock gas O2. In the first case, rij is the number 
of simulated particles of the respective type in the con- 
sidered cell divided by the cell volume, whereas in the 
second case, nj = p/kT, with p the gas pressure and T 
the room temperature. We do not calculate the mean 
free times tj for all the collisions listed in tabic HI In- 
stead, we set Tj = Ate, 2 for Coulomb collisions (1) - (3) 
and Tj = Aic,i for the remaining collisions (4) - (19). 
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TABLE I: Collisions included in our model 

elastic scattering 

(1) e + e ^ e + e 

(2) 0~ +0" ^ +0- 

(3) 0+ + 0+ ^ 0+ + Ot 

(4) e + Oa ^e + Oz 

(5) +O2 ^ O" +O2 

(6) 0+ + O2 ^ O2 + 0+ 
electron energy loss scattering 

(7) e + 02 ^e + 02(/^ = 1,...,4) 

(8) e + 02 ^e + 02(Ryd) 

(9) e + O2 ^ e + 0(3P) + 0(3P) 

(10) e + O2 e + 0(3P) + 0(1D) 

(11) e + 02 ^e + 02(aiAg) 

(12) e + 02 ^ e + 02(b^Eg) 
electron & ion production & loss 

(13) e + Oj ^ + 

(14) O- + 0+ ^ O + O2 

(15) e + 02 ^ + 0" 

(16) O" + O2 ^ O + O2 + e 

(17) O- +02(a^Ag) -» O3 +e 

(18) e + O2 ^ 2e + Oj 

(19) e + O^ ^ + 2e 



charge exchange 

vibrational excitation 
Rydberg excitation 
dissociation (6.4 eV) 
dissociation (8.6 eV) 
meta-stable excitation 
meta-stable excitation 

dissociative recombination 
neutralization 
dissociative attachment 
direct detachment 
associative detachment 
impact ionization 
impact detachment 



In order to verify that our code treats collisions cor- 
rectly, we performed various test runs and compared the 
results either with analytical results 45| or with results 
obtained from the BITl code [s^. The latter was also 
used to check discrepancies in the cross section data. 



III. SCATTERING AND REACTION CHANNELS 

Collisions strongly influence the particle concentration 
and the energy balance in the discharge. For oxygen, an 
overwhelming number of elastic, inelastic, and reactive 
collisions is possible. Restricting, however, the descrip- 
tion to a three species plasma (e, O^, and Oj) many 
processes can be ignored and others can be treated ap- 
proximately. In tableUwe show the collisions defining our 
three species model for an oxygen discharge, and classify 
them into three groups: elastic scattering, electron en- 
ergy loss scattering, and electron and ion production and 
loss reactions. The respective cross sections for these 
processes are shown in figures |3] - [S] as a function of the 
relative energy of the two particles in the entrance chan- 
nel of the respective process. 

Our collection of cross sections is semi-empirical, com- 
bining measured data, which is usually available only in 
a finite energy range, with simple models for the low- 
energy asymptotic, which in most cases is not very well 
known from experiments. In general, the high-energy 



asymptotic has to be also determined from models, but 
it is less critical because the distribution functions usu- 
ally decay sufficiently fast at high energies. If not stated 
otherwise, we extrapolated therefore the values of the 
cross sections for the largest energies shown in the plots 
to all energies above it. Some of the cross sections signif- 
icantly deviate from the ones previously used [s^. Our 
simulations indicate, however, that the modifications are 
essential for obtaining bulk densities in accordance with 
3. 



experiments 
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A. Elastic scattering 



Elastic collisions (1) - (6) are particle number con- 
serving. Thus, when a collision takes place, only the 
post-collision velocities of the scattering partners have 
to be determined whereas the list of simulated particles 
remains unaltered. 

The first three scattering processes (1) - (3) are intra- 
species Coulomb collisions. They are not very important 
for the experiments 27, 13] we analyze in this series of 
papers. For other parameter regimes, however, the local 
ion density in an electro-negative gas discharge can be 
rather high ^23*1 and it cannot be ruled out that ion-ion 
Coulomb collisions affect, for instance, the width of these 
ion density peaks. Since the computational burden is 
moderate, we kept intra-species Coulomb scattering in 
all our simulations. If not stated otherwise, inter-species 
Coulomb collisions were however neglected. 

For Coulomb collisions, we used a binary collision 
model 44] with a uniformly distributed azimuth angle 
0, and a Gaussian distribution for tanx/2, where x is 
the scattering angle. The second moment of the Gaus- 



sian distribution is given by 44 ] 



(tan' 



2 



q^n In F 



tr 



(2) 



with tc,2 the collision time for Coulomb scattering, InF 
the Coulomb logarithm, u the magnitude of the relative 
velocity in the center-of-mass frame, and n, g, m, and 
eo the local density, the charge, the reduced mass of the 
charged species under consideration, and the dielectric 
constant of the vacuum, respectively. 

Figure [3] shows the cross sections for elastic scattering 
of electrons (4) and 0~ ions (5) on O2 molecules together 
with the cross section for charge exchange scattering (6) 
of O^ ions on O2 molecules. Since little is known about 
the angle dependence of these processes, we assume the 
collisions to be isotropic, with a uniform distribution of 
the azimuth angle 4> and a uniform distribution of cos x, 
where x is again the scattering angle. 

The cross section for elastic scattering of electrons on 
O2 [il] is the same as in Ref. [sl]. The other two cross 
sections are different. For charge exchange scattering, we 
constructed the cross section as follows. Below 0.251 eV 
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FIG. 3: Cross sections for elastic scattering of electrons (4), 
0~ ions (5), and ions (6) on oxygen molecules. The 
dashed line and the dotted line indicate, respectively, the cross 
section for elastic (0~ , O2) scattering and the cross section for 
(0^,02) charge exchange scattering used in Ref. [s^. 



and above 8.5 eV, we used empirical data for momentum 
transfer scattering 4^, H^l , together with the expression 



(3) 



where cr™ and <7cx denote, respectively, the momentum 
and charge exchange cross section. For energies in be- 
tween, we employed a linear interpolation. With this 
cross section, we could reproduce experimentally mea- 
sured O2 velocity distribution functions very well 1431 . 
With the charge exchange cross section given in Ref. [36[ 
(dotted line in figure [3]), on the other hand, we could not 
obtain the correct distribution functions - neither with 
our PIC-MCC code nor with the BITl code |385, which 
we used in addition to exclude possible mistakes in our 
collision routines. 

The cross section for elastic scattering of on O2 
originates from a model which assumes polarization scat- 
tering between the two, induced by a central potential 
V{r) — ~C/r" with n — 3. It is given by 



<^e{E) = [- 



C 



,2/7^ 



E[eV] ■ 



(4) 



with C = 3.77-10-24 and do = 1-2 -IQ-^^ Here, 
we added a constant shift (Tq to match the cross section 
at high energies as given in Ref. 36]. 



B. Electron energy loss scattering 

Electron energy loss occurs due to inelastic collisions, 
(7) - (12), in which the oxygen molecule is either pro- 
moted into an excited state or dissociated into neutral 
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excitation of O2 molecules (7). 



fragments. In the three species model, where the spatio- 
temporal evolution of the neutral fragments and the ex- 
cited oxygen molecules is not traced, we can treat this 
group of processes in the spirit of "test particle collisions" 
where the electron is the "test particle", suffering mo- 
mentum as well as energy transfer, and the O2 molecule 
is the "field particle" (with internal degrees of freedom). 
Drawing the field particles from a Maxwell distribution 
characterized by the gas temperature T and the gas den- 
sity n = kT/p, the same binary collision model can be 
used as for elastic scattering [4J]. The only difference 
is that now the magnitude of the post-collision relative 
velocity in the center-of-mass system is given by 



2SE 

mafj 



(5) 



where u is the magnitude of the pre-collision relative ve- 
locity in the center-of-mass frame, rriap is the reduced 
mass of the scattering partners, which, in general, have 
not the same mass, and 5E is the excitation energy of the 
process. Due to lack of angle-resolved scattering cross 
sections, we again assume the collisions to be isotropic. 

The most important electron energy losses are due to 
vibrational and electronic excitation and impact dissoci- 
ation of the oxygen molecule. The cross sections for these 
processes [il] are the same as the ones used in Ref. [s^. 
For completeness, they are plotted in figures S] and [H 
We consider vibrational excitations with v = 1,2,3, and 
4, excitation of the O2 Rydberg state, PP and PD dis- 
sociation, and the excitation of the meta-stable states 
02(a^Ag) and 02(b^Sg). In contrast to Ref. [11], we ig- 
nore DD dissociation, because, in the energy range of in- 
terest, its cross section is one order of magnitude smaller 
than for PP and PD dissociation, and rotational excita- 
tion which is almost elastic (small energy transfer) and 
thus dominated by the much more efficient elastic scat- 
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FIG. 5: Cross sections for Rydberg excitation (8), 6.4 eV (9) 
and 8.6 eV (10) dissociative excitation, A excitation (11), and 
E excitation (12) of O2 molecules. 



tering of electrons on O2. 
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FIG. 6: Cross sections for dissociative recombination (13), 
neutralization (14), dissociative attachment (15), direct de- 
tachment (16), associative detachment (17), impact ionization 
(18), and impact detachment (19). The dashed line indicates 
the experimentally determined high-energy asymptotic of the 
neutralization cross section and the dotted line is the cross 
section for neutralization used in Ref. [S^. 



C. Electron and ion production and loss reactions 

Electron and ion production and loss reactions (13) - 
(19) are inelastic collisions which change the number of 
charged particles. Therefore, the list of simulated parti- 
cles has to be up-dated. All coUisions obey conservation 
laws and are again assumed to be isotropic. The details 
of the modeling depend on the process. 

The two recombination channels (13) and (14) result 
simply in the annihilation of the two oppositely charged 
particles which participate in the process. Binary colli- 
sions with one charged particle in the exit channel, such 
as dissociative attachment (15) and associative detach- 
ment (17), are treated within the modified binary colli- 
sion model described in the previous subsection. Impact 
ionization (18) and impact detachment (19) are mod- 
elled as follows: First, an inelastic binary collision is per- 
formed, in which the parent electron looses the ionization 
(detachment) energy. Then, the post-collision O2 (0~) 
particle is split into an electron and a (O) particle. 
Finally, an elastic binary collision is applied to distribute 
energy among the two charged particles. Modeling three 
particle processes as a sequence of two binary collisions 
with a particle splitting in between guarantees energy 
and momentum conservation which is critical for the sta- 
bility of simulations 45|. Although direct detachment 



(16) could be modelled in the same spirit, we adopted a 
simpler approach. From experiments [gO] we know that 
the energy of the released electron is approximately one- 
tenth of the energy of the primary 0~ ion. To obtain the 
energy distribution of the ejected electron we hardwired 
therefore this ratio in the collision routine for direct de- 
tachment. 



The cross sections for electron and ion production and 
loss reactions are shown in figure [HI For impact ion- 
ization [4^, dissociative attachment [521] . and impact 



detachment [53|, l54|, |55| we used empirical cross sec- 



tions throughout, putting, in the case of impact detach- 
ment a higher confidence level to more recent data [55| . 
whereas for dissociative recombination, ion-ion neutral- 
ization, and detachment on neutrals we combined empir- 
ical data with analytical models. 

First, we discuss dissociative recombination (13). For 
this process we used an effective cross section. 



ad,iE) ^ B ■ aoiE) 



(6) 



where ctq is the cross section for non-resonant dissociative 
recombination (63l | and = 3 is a scaling factor, which is 
used to fit experimental data [1^. A comparison of the 
effective cross section with a calculated cross section [63 | 
shows that B can be also considered as a simple device to 
approximately take the effect of resonant dissociative re- 
combination channels into account which originate from 
vibrationally excited states of . 

For ion- ion neutralization (14), we constructed a cross 



section from a two-channel Landau-Zener model 6l|, 
with one free parameter, which we adjusted to obtain the 
correct high-energy asymptotic of the cross section [st} . 
Ion-ion neutralization occurs because the adiabatic en- 
ergy of the (0^,0^) configuration decreases when the 
ions approach each other. At a certain distance Rx, the 
energy of the (0^,0^) configuration falls below the en- 
ergy of the (O, O2) configuration [65[. The Landau-Zener 
theory estimates the probability for switching from one 
configuration to the other at the distance and leads 



to a cross section, 



ariE) = 47ri?^(l 



RxE 



(7) 



which approaches a constant C = 47r_R^ for ^ 1- 

Fitting C and thus Rx to empirical data at high ener- 
gies [57|, we obtain 



a.(i?) = 0.8.(l + |^).10-^.cn^^ 



(8) 



This cross section deviates dramatically from the one 
used in Ref. (dotted Hue in figure 0- 

The most severe modification we made was for detach- 
ment of 0~ on neutrals. It takes place through direct 
detachment (16), resulting in an atomic oxygen, an oxy- 
gen molecule and an electron, and through associative 
detachment (17), which leads to an O3 molecule and an 
electron. The latter is rather surprising because there 
is no evidence for it in beam exp eriments (where only 
direct detachment is observedJ[5gj59| ) . Yet, experimen- 
tal studies of O2 discharg es 1271. 1681 fas well as general 
theoretical considerations |l6l. 1661. l67j) strongly suggest 
that associative detachment is possible in an oxygen dis- 
charge because of the presence of meta-stable ©2(0^ A^). 
In contrast to direct detachment, which has a threshold 
around 1.3 [11] (see figure [6|), associative detachment 
has usually no threshold. Thus, it may be a rather im- 
portant loss channel for cold 0~ ions. 

Since we could not find an empirical cross section for 
associative detachment, we employed a model, which de- 
scribes the detachment (electron loss) as the "inverse" of 
a classical Langevin-type electron capture into an attrac- 
tive auto-detaching state of O;^. Assuming the polariz- 
ability for 02(0^ Ag) to be the same as for O2, the cross 
section is then given by [65| 



a^^{E) = 5.96 • 
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-16 



(9) 



As can be seen in figure [6l associate detachment is al- 
ready the dominant detachment process for energies be- 
low ~ 6 eV . 

From the cross section alone, of course, we cannot de- 
termine the probability Pad for associative detachment. 
We also need the density of 02(a^Ag), which is unknown 
in the three species model. However, it should be of the 
order of the O2 density. Therefore, we write tt-a = C-uq^, 
with C < 1, and obtain Pad = u ■ cFad • (C* • nQ^) ■ t^^i, 
where C is a fit parameter which can be adjusted to ex- 
periments. 



IV. COMPARISON WITH EXPERIMENT 

The plasma model described in the previous sections 
contains three parameters: ctq, B, and C. The first two 
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FIG. 7; Voltage dependence of the quasi-stationary, cycle- 
averaged electron and ion densities at x = L/2 (bulk den- 
sities) for a discharge with L = 2.5 cm, p = 13.8 Pa, and 
frj ~ 13.6 MHz. Filled and open symbols denote, respec- 
tively, results of the simulation and experimentally measured 
densities [13 . Solid and dashed lines guide the eye. 



we fixed by a direct comparison of the model cross section 
with the measured cross section. The parameter C, in 
contrast, cannot be determined in this way, because it is 
not linked to a cross section. It denotes instead the frac- 
tion of O2 molecules in the meta-stable 02(a^Ag) state 
which in turn depends on the particular set-up of the 
discharge. We consider therefore C as a free parameter 
which can be adjusted to the discharge to be modelled. 

To validate our model we focus now on the oxygen dis- 
charge described in Ref. 27]. In order to determine C, we 
simulated the discharge for p — 13.8 Pa, Urf — 250 V, 
and /,./ = 13.6 MHz, and tuned C to reproduce the 
quasi-stationary, cycle-averaged negative and positive ion 
densities at x = L/2 (bulk densities); the bulk electron 
density matches then also because of quasi-neutrality in 
the bulk of the discharge. We obtained C w 1/6, imply- 
ing that roughly one sixth of the O2 molecules is in the 
meta-stable state. 

In figure [7] we plot the bulk ion and electron densities 
for p — 13.8 Pa over a wide voltage range, using however 
for all voltages C ~ 1/6, the value determined for Urf = 
250 V. As can be seen, the agreement between simulation 
and experimental data (27j is rather good, indicating that 
our model captures the essential processes in the bulk of 
an oxygen discharge. 

The precise value of C should not be taken too serious 
because it is based on a rather crude model for associa- 
tive detachment. More important is that without this 
process (C = 0), the simulation could not reproduce the 
measured densities. This can be seen in figure [H where 
we plot the quasi-stationary, cycle-averaged density pro- 
files without (thin lines) and with associative detachment 
(thick lines) taken into account. (For these runs, we also 
included inter-species Coulomb collisions but the profiles 
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FIG. 8: Quasi-stationary, cycle- averaged electron and ion 
density profiles for U = 150 V without (thin lines) and with 
(thick lines) associative detachment taken into account; p, L, 
and frf are the same as in figure [T] Open circles and open 
squares denote, respectively, the measured electron and 0~ 
density profiles [23|. Note, the runs producing the shown data 
included inter-species Coulomb scattering. 



without it are basically the same.) Neglecting associative 
detachment, the bulk ion densities turned out to be al- 
most one order of magnitude too high. Remarkably, with- 
out associative detachment, the density profiles have the 
parabolic shape expected from ambipolar drift-diffusion 

models [13, Ell, El, til, Q H Ei ■ 

To understand why associative detachment is so cru- 
cial for modeling the discharge of Ref. [23|, we plot in 
figure[9]for p = 13.8 Pa and Urf = 150 V the (marginal) 
velocity distribution function for negative ions, f-{E,x) 
with E = sign{vx)vl/2M^. Clearly, the majority of 
ions is cold; the distribution function spreads out only 
in the region where the bulk plasma merges with the 
peripheral plasma at x ~ 0.8 cm. Thus, loss processes 
whose cross sections are large for small energies (see fig- 
ure [6|), that is, ion- ion neutralization and associative de- 
tachment, will be very efficient. Direct detachment, on 
the other hand, whose cross section has a threshold, will 
be suppressed. At low energies, both ion-ion neutral- 
ization and associative detachment have cross sections 
which increase with decreasing energy. Which process 
dominates depends therefore also on the collision prob- 
ability P. For ion-ion neutralization P is proportional 




FIG. 9: Quasi-stationary, cycle-averaged negative ion velocity 
distribution function f-{E,x) with E = sign{vx)Vx/2M- for 
Urf = 150 V; p, L, and frf as in figure [3 



that associative detachment is the main loss process for 
0~ ions. First, we conclude from figure flOl which shows 
the attachment cross section (in arbitrary units) together 
with bulk electron velocity distribution functions for dif- 
ferent voltages, that the energy range between 4 eV and 
15 eV will be the most important one for the attachment 
process; for energies below 4 eV the production is zero 
because the cross section is zero while for energies above 
15 eV the number of electrons available for attachment 
is too low. We also see that in the voltage range of inter- 
est, the electron velocity distribution function does not 
change much in this energy range. Thus, the production 
process is almost voltage-independent and the observed 
near constancy of the 0~ density has to be taken as 
a signature of the loss process. As already mentioned, 
direct detachment can be ruled out. Ion- ion neutraliza- 
tion, on the other hand, would lead for E <^ lA eV, that 
is, in the relevant energy range, to an energy-resolved 
rate coefficient Kr{E) ^ ypE ■ ar(E) ^ 1/VE which 
strongly increases with decreasing energy. Thus, even 
small, voltage-induced changes of the energy of negative 
ions would suffice to lead to noticeable modifications of 
the loss rate. What remains is associative detachment. 
The energy-resolved rate coefficient for this process is 



to the O, density while for associative detachment it is 
proportional to the 02(a^Ag) density. Since for the dis- 
charge of Ref. [27|, the density is much smaller than K^^{E) ^ \fE ■ cr^j{E) ^ const. At the same time, the 

distribution function for ions does not change much 
between 75 V to 300 V. Thus, the loss rate due to as- 
sociative detachment is nearly independent of voltage. 
Together with the voltage independence of dissociative 
attachment, this leads to the observed constancy of the 
0~ density. 

Our simulation reproduces the absolute values of bulk 
densities fairly well. The (axial) density profiles of the 



the 02(a^Ag) density, which we estimated to be 1/6 of 
the O2 density, associative detachment has to be the main 
loss process for negative ions in this experiment. 

Additional support for this claim comes from the weak 
voltage dependence of the 0~ density between 75 V and 
300 V indicating that both 0~ production and loss 
are in this voltage range rather insensitive to voltage 
changes. We will now show that this strongly suggests 
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FIG. 10: Quasi-stationary, cycle-averaged electron velocity 
distribution function fe{E,x) with E = w^/2me and x — L/2 
for Urf = 75 V,150 V,280 V, and 300 V; p, L, and frf as 
in figure [T] The dotted line indicates the cross section for 
dissociative attachment (15) in arbitrary units. 

simulation decay however too fast compared to the ex- 
perimental ones (open symbols in figure[5]). We attribute 
this to the fact that the simulations assume a constant 
02(0"'^ Ag) density, whereas, in reality, there should be 
a density profile, because of the interplay of volume 
and surface loss and generation processes for 02(a^Ag) 
molecules. In particular surface losses, that is, the decay 
of ©2(0^ Ag) molecules when they hit the boundary of 
the discharge should play an important role because they 
lead to a depletion of □2(0"'^ Ag) molecules in the vicinity 
of the electrodes and thus to a reduction of associative 
detachment. A detailed investigation of the "drag" the 
©2(0^ Ag) density is expected to exert on the 0~ density 
is however beyond the scope of the present paper. 



V. CONCLUSIONS 

We presented a planar, one-dimensional PIC-MCC 
model for a capacitively coupled rf discharge in oxygen 
which is capable to quantitatively describe experiments. 
Its main features are: (i) Only electrons, 0~ ions, and 
O2 ions are treated kinetically. (ii) Neutral particles are 
only incorporated in as far as they affect the particle and 
energy balance of simulated charged particles, (iii) A 
direct simulation Monte Carlo model for collisions has 
been used, with two groups of collisions, which differ in 
the collision times. 

The elementary processes of our three species model 
are shown in table ID We took the processes with the 
largest cross sections into account. When possible, we 
implemented empirical cross sections. In some cases, 
however, we combined them with analytical model cross 
sections to obtain the correct low-energy asymptotic. 
In particular, we used a polarization-type scattering 



model for elastic (0~,02) scattering, a Landau-Zener- 
type model for (0~,0j) neutralization, and an inverse 
classical, Langevin-type capturing model for 02(a^Ag)- 
induced associative detachment. To describe associative 
detachment due to ©2(0^ Ag) without explicitly calcu- 
lating the ©2(0^ Ag) density, we furthermore introduced 
a parameter C, which should be fitted to the particular 
oxygen discharge under consideration and could be inter- 
preted as the fraction of O2 molecules in the meta-stable 
02(0"^ Ag) state. 

As a first application of our model, we simulated the 
discharge of Ref. [13] ■ After we adjusted the parameter 
C to reproduce the bulk ion densities for p = 13.8 Pa 
and Urf — 250 V, we also obtained the correct bulk den- 
sities for other voltages. The parameter C « 1/6, that is, 
approximately one sixth of the oxygen molecules are in 
the meta-stable state. We also pointed out that the weak 
voltage dependence of the bulk densities for fixed pres- 
sure indicates that 0~ losses due to associative detach- 
ment dominate the ones due to ion-ion neutralization. 

Although we could reasonably well describe electron 
and ion densities in the bulk, the (axial) ion density pro- 
files of the simulation are too narrow compared to the 
experimental ones. Most probably this is because we as- 
sumed a constant 02(0"'^ Ag) density. In reality, however, 
02(a^Ag) molecules decay when they hit the boundary 
of the discharge. The 02(0^ Ag) density should there- 
fore decrease in the vicinity of the electrodes and with it 
associative detachment. To take this effect into account 
requires however a model which not only treats electrons, 
0~ ions, and ions kinetically, but also 02(0^ Ag) 
molecules. 

Our concern in the present paper, which is the first 
in a series of three, was the calibration and validation 
of a three species model for capacitively coupled rf dis- 
charges in oxygen. For that purpose, we used published 
experimental data for the bulk of the discharge [27| . 
The following paper js^l describes the results of an 
experimental investigation of the sheath region of such 
a discharge, focusing, among others, on the spatial and 
temporal evolution of the 844 nm double emission layer 
in front of the powered electrode. In the third paper [43(, 
finally, we will use the three species model to simulate 
the sheath region. Besides a comparison of simulated 
with measured distribution functions for positive ions, 
we will also present a model for the 844 nm double 
emission layer. The three species model, including 
its approximate treatment of associative detachment, 
is sufficient for that purpose, because negative ions 
are negligible in the region from which the emission 
originates. Negative ions are only needed for the overall 
charge balance of the discharge. Their density profile 
per se is not crucial for the explanation of the 844 nm 
double emission layer as long as it features a vanishingly 
small density in front of the electrode. This is already 
accomplished by the three species model. 
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